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ABSTRACT 

Semi-regular (SR) variables are not a homogeneous class and their variabil¬ 
ity is often explained due to pulsations and/or binarity. This study focuses on 
IRAS 19135+3937, an SRd variable with an infra-red excess indicative of a dusty 
disc. A time-series of high-resolution spectra, UBV photometry as well as a very ac¬ 
curate light curve obtained by the Kepler satellite, allowed us to study the object 
in unprecedented detail. We discovered it to be a binary with a period of 127 days. 
The primary has a low surface gravity and an atmosphere depleted in refractory el¬ 
ements. This combination of properties unambiguously places IRAS 19135+3937 in 
the subclass of post-Asymptotic Giant Branch stars with dusty discs. 

We show that the light variations in this object can not be due to pulsations, 
but are likely caused by the obscuration of the primary by the circumbinary disc 
during orbital motion. Furthermore, we argue that the double-peaked Fe emission lines 
provide evidence for the existence of a gaseous circumbinary Keplerian disc inside the 
dusty disc. A secondary set of absorption lines has been detected near light minimum, 
which we attribute to the reflected spectrum of the primary on the disc wall, which 
segregates due to the different Doppler shift. This corroborates the recent finding that 
reflection in the optical by this type of discs is very efficient. The system also shows a 
variable Ha profile indicating a collimated outflow originating around the companion. 
IRAS 19135+3937 thus encompasses all the major emergent trends about evolved disc 
systems, that will eventually help to place these objects in the evolutionary context. 

Key words: circumstellar matter - stars: AGB and post-AGB - binaries: spectro¬ 
scopic - stars : variables: general 


1 INTRODUCTION 

There are four types of semi-regular (SR) variables, each 
designated alphabetically from a to d. IRAS 19135+3937 
belongs to the last category, "d". SRd variables represent 
a poorly understood group of warm (spectral types F-K) 
giants and supergiants. Their light curves exhibit variations 
on the time-scales 30-1100 days, that are only quasi-regular 
and lack characteristic features of the radial pulsators with 
similar periods, such as Cepheids and RV Tau stars. There 
is, however, one intriguing feature in the behaviour of some 
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RV Tau stars that is reminiscent of the SR variability. Nor¬ 
mally RV Tau stars pulsate with periods between 30-150 
days, with two alternating minima per cycle, shallow and 
deep. An ’RV Tau b’ subgroup, however, shows an additional 
variability in the mean magnitude, amplitude, or relative 
strengths of the minima. This additional ’secondary period’ 
variability occurs on the time scale that overlaps with vari¬ 
ability in the SRd-s with long periods. This phenomenon is 
referred to as ’long secondary periods’ (L SP), and the ori¬ 
gin of LSPs is still a matter of debate fe.g.. lKiss et al.ll2007l : 
I Wood fe et al .1 1 1 9991 ) . The two major hy potheses discussed 
in the litera t ure are variable obs curation (ILlovd E van J 1 9741 : 
iFokinl [l994l : IPollard et al ]ll997|) and an in terpla y between 
various pulsation modes dBuchler fe Kovacslll987h . 
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Modelling of the better studied Galactic systems is hin¬ 
dered by the lack of direct distance measurements, hence 
luminosities. In terms of kinematics, chemical composition, 
and emission line strength there appears to be no sy stem ¬ 
atic d ifference between SRd-s and RV Tau stars JWahlgrenl 
GUI). Chemical composition studies indicate that both 
groups are heterogeneous and include halo, as well as disc 
objects rtGiridhar et aLlfeoOOUBritavskiv et al Jl2()T3 l. Based 
on the rarity, presence of the circumstellar matter, and 
the luminosities deduced for the Magellanic Clouds ob¬ 
jects, most of these objects must represent late evolutionary 
stages of low mass stars, being either post -red g iant or post- 
asymptotic giant branch (pAGB) stars dWallersteinl 12002 : 

I Van Winckel! 120031 : iKamath et al.l [20141 ). Since there is no 
evidence that SRd-s and RV Tau stars would represent dif¬ 
ferent stellar populations or different evolutionary stages, 
the reason for the semi-regular behaviour must be sought in 
the properties of individual systems. 

Binarity could be such a property, but how exactly 
would it cause the observed light variations? Mutual eclipses 
of the companions can be ruled out, as the fraction of sys¬ 
tems with an edge-on orientation or where both compan¬ 
ions are evolved giants, must be negligible. Furthermore, 
the semi-regular character requires that the eclipsing body 
should be variable. It has long been known that many RV 
Tau stars have near-infrared excesses, indicating presence o f 
hot dust up to the sublimation temperatures llEvan 
The modelling of the spectral energy distribution (SED) 

( Gielen e t al ] 200! ) , as well as interfero metric measurements 
( Deroo et al.l 2007 ; iHillen et al.l I2013T ) have demonstrated 
that the dust resides in a disc. The discs have typical in¬ 
ner radii of only a few astronomical units (~10 stellar radii, 
R*) and inner walls of substantial scale-height. 

At the same time, there is a growing evidence 
from the radial velocity (RV) studies that evolved ob¬ 
jects with discs are all (binaries llVan Winckel et al.1119991 : 
IVan Winckel fc et alJ 12009 1. Based on the deduced orbits 
for the visible giant primaries, the discs must be circumbi- 
nary. The variability can then be naturally explained by the 
obscuration of the primary star by the inner rim of the disc, 
as the line of sight to the star pr obes different heights abov e 
the disc plane along its orbit llWaelkens fe Watersl 1 1993lf . 
The cycle-to-cycle irregularities could be caused by inhomo¬ 
geneity or precession of the disc. 

To test this theory, one has to find a correlation between 
the semi-regular behaviour, the shape of the SED, and bi¬ 
narity. In reality, it is not always possible to prove binarity 
based on RV variations alone, for example in case of small 
inclination angle or strong pulsations. In 2009 we started to 
monitor a number of RV Tau, SR, W Vir, and chemically pe- 
culiar pAGB stars with the echelle spectrogra ph HERMES 
llVan Winckel et al.l l2010l : iGorlova et al.l 1201 il l . Besides de¬ 
tecting RV variations consistent with binarity and obtain¬ 
ing orbital solutions for some systems, in many of them 
we also discovered specific phase-dependent Ha profiles 

l|Gorlova et al.ll2012l: IVan Winckel et akll2012l: IGorlova et al.l 

120131! . Normally these profiles are P Cyg-like or have a 
double-peak emission, but only during the superior conjunc¬ 
tion of the pAGB prima ry they were found to develop strong 
blue-shifted absorption. IWitt et al.l ll2009l l. based on the ob¬ 
servations of a similar phenomenon in the central star of the 
Red Rectangle (RR) nebula, explained the Ha line-profile 


variability by a model in which a jet is powered by accre¬ 
tion from the giant primary to the invisible (likely a main- 
sequence, MS) companion. IRAS 19135+3937, originally in¬ 
cluded in our sample due to the disc-like infra-red (IR) ex¬ 
cess, turned out to be one of such objects. Besides Ha, it 
caught our eye also because of the rather smooth light-curve 
with a large amplitude, which is not typical for long-period 
pulsators. IRAS 19135+3937 thus presented a good case for 
testing binary theory for semi-regular variables. 

The variability of IRAS 19135+3937 has been first dis¬ 
covered by amateur astronomers (ISallman fe Droegell2004ll . 
who determined a period P=125.4 d and an amplitude 
AV = 0.9 mag. Based on these observations the star was 
included in the General Catalogue of Variable S tars as an 
SRd variable V677 Lyr llKazarovets et all 1 201 .'j |. The star 
was also observed in the course of the All Sky Automated 
Survey (ASAS), where it was classified as a “QPER” type, 
which i s a sem i-r egula r variable with a stable period (128.8 
d. lPigulski et al.ll2009l f. The star is relatively faint (V ~ 11 
mag) and had not been discussed i n the literature when we 
started its observations. Recently, iRao fe Giridharj (|2014l f 
performed an abundance analysis and found it to be a mod¬ 
erately metal-poor star ([Fe/H]=—1 dex) with some pecu¬ 
liarities_ 

In lGorlova et al.1 l|2012l l we presented the first RV curve 
of IRAS 19135+3937 that revealed its binarity. We also 
showed a number of peculiar spectral lines including Ha, 
and pointed to the similarity with BD+46°442, the first ob¬ 
ject discovered in ou r survey to displa y this type of Ha line 
profile variability IlGorlova et al.ll2012h . Since then we have 
tripled the number of observed spectra of IRAS 19135+3937 
and followed up with multi-band photometry. Here, we pro¬ 
vide a comprehensive analysis of these data in order to un¬ 
derstand the cause of the semi-regular variability in this 
object. We deduce strong observational constraints on the 
geometry of the system and study the ongoing interaction 
processes. 


2 NEW OBSERVATIONS AND DATA 
REDUCTION 

2.1 Photometry 

We observed IRAS 19135+3937 with a photoelectric UBV- 
photometer dLyutvl| [T97ll) attached to the 60-cm telescope 
Zeiss-1 (located in Crimea) of the Sternberg Astronomical 
Institute (SAI Russia). The diameter of the photometer di¬ 
aphragm was set to 27 arcsec. Seventy four measurements 
have been made between 2012 - 2013, which overlaps with 
the epoch of our spectroscopic observations. 

To obtain colour terms for the transformation of the in¬ 
strumental ubv system to the standard Jo hnson’s UB V one, 
we observed standard stars in NGC 6633 l)Mermilliodll 198(f ). 
By solving equations AV = Av — kiA(b — v), A(B — V) — 
k 2 A(b — v), and A (U — B) = k 3 A(u — b), we obtained the fol¬ 
lowing colour terms: k\ = 0.082 ± 0.010, k 2 = 0.980 ± 0.010, 
fc 3 = 1.081 ± 0.015. We used HD 179909 as a comparison 
star, with the following magnitudes adopted from Mermil- 
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liod (19910: V = 8.285 ± 0.015, B-V = -0.070 ± 0.010, 
U—B = —0.330±0.010. The typical accuracy of our absolute 
photometry is of the order of 0.01 mag. 

The measurements are given in Table EH The obser¬ 
vations revealed sinusoidal-like variations in brightness and 
colour, with a period consistent with earlier reports. 


2.2 Spectroscopy 

We collected 61 spectra of IRAS 19135+3937 with the HER¬ 
MES fibre echelle spectrograph m ounted o n the 1.2 m tele¬ 
scope Mercator on La Palma ll Raskin et al.l fzoTll i . The ob¬ 
servations were carried out in the years 2009 - 2013, with 
an average cadence of one observation per two weeks. The 
exposure time varied with the brightness of the star and 
the weather conditions, being on average 1400 s. A typical 
signal-to-noise ratio (S/N) in the central orders (around Ha) 
was ~35, which was needed for obtaining a cross-correlation 
function suitable for the RV determination (Section 15.11) . 
All spectra in our survey have been collected using the high- 
resolution fibre configuration, that provides R ~ 80, 000 and 
a wavelength coverage between 3800 and 9000 A, and re¬ 
duced with a dedicated PYTHON-based pipeline Hermes- 
DRS. 

The pipeline provides science graded output. A quick- 
look of the extracted spectrum is available during the ac¬ 
tual observations, and the full reduction is automatically 
performed after sunrise. The pipeline first averages in 2D 
biases, arcs, and flat fields. The cross-order profile with the 
two slices is modelled on a daily basis and used at extraction, 
during which the cosmic clipping is also performed. We do 
not use optimal extraction. The read-out noise is typically 
3-4 electrons. Sky subtraction is not performed, because all 
spatial information is lost due to the scrambling properties 
of the optical fibre. The contamination by the interstellar 
emission, however, is not a concern for our evolved stars, 
that are located away from star-forming regions. 

The stability of the zero-point of our arc-based wave¬ 
length calibration is monitored with the help of the Inter¬ 
national Astronomical Union RV standards. Over 5 years 
of HERMES operation we obtain a standard deviation of 
0.080 km/s based on 2329 measurements of different radial- 
velocity standards, measured as a spread over the mean of 
every standard. The shifts are mainly cased by t he pressure 
variations during the night jRaskin et al.l [2011 1. More de¬ 
tails on the instrument and dat a reduction can be fo und at 
the HERMES webpagc0 and in lGorlova et al.l (l2012l b 


3 BASIC PROPERTIES 

3.1 Photospheric parameters and chemical 
composition 

For the chemical abundance study we used an average of 
three well exposed consecutive spectra obtained on 29 June 
2009, with maximum S/N~150 near 6000 A in the combined 


1 VizieR on-line catalogue 11/168 “Homogeneous Means in the 
UBV System”. 

2 http://www.mercator.iac.es/instruments/hermes/ 


spectrum. These observations were carried out near maxi¬ 
mum light, and metallic lines were nearly symmetric and 
easy to measure. As will be shown in Sec. 15.11 this phase 
corresponds to the superior conjunction of the pAGB pri¬ 
mary, when its obscuration by the circumstellar disc and the 
contribution from a putative companion should be minimal. 

We followed the same procedure for the physical pa¬ 
rameters and abundance determination as for a n other sim¬ 
ilar object BD +46°442, as described in iGorloval (1201 il l and 
1 Gorlova et al.l ( 2012? !. Briefly, to get the initial estimate of 
the effective temperature (T e ff) and surface gravity (logg), 
we matched the observed profiles of H<5, Hy, and H/3 (wings 
only) with the prec omputed state of the art model pro- 
files of Coelho et al. (12 00511 . and Paschen 14 with those of 
iMunari fe Castelli! ( 200(1 . The two distinct sets of model 
spectra had to be employed in order to cover both Balmer 
and Paschen series. After visual inspection^ we selected the 
following best combinations of T/fj/log 6000/1.0+0.5, 


6250/2.0, and 5750/0.5. 

At the next stage of iteration, we examined the cor¬ 
relation of the Fe abundance obtained from the Fel and 
Fell lines with their equivalent widths (EWs). The abun¬ 
dances (for Fe, as well as other elements) were computed 
from the EWs using MOOGIO, which is the latest ver¬ 
sion of the L TE abundance determination code by C. Sne- 
den dSnedenl Il973ll . and the ATLAS9 model atmospheres 
of R. Kurucz in the updated version of ICastelli fe Kuruczl 
lj2003i l. We started with models with the solar metallicity, 
but re-computed the final abundances with [M/H] = —1.0 
to better match the deduced overall metal deficiency of 
IRAS 19135+3937. The EWs in the observed spectrum 
were measured using our PYTHON-based program, which 
includes a module to disentangle close line pairs. The lat¬ 
ter is important for IRAS 19135+3937, because the lines 
are relatively broad (FWHM = 12.5 km s -1 in the consid¬ 
ered phase). To maximize the number of measured lines, 
we also augmented our pre vious atomic li ne list ( which was 
an up-dated version from iKovtvukh fe Andrievsky! fl999h . 
with two others: a list previou sly used by the Leuven group 
llVan Winckel fc Revniersll2000h. which is largely b ased on 
the solar list by F. Thevenin jTheveninlll989l . 1 199dh , and a 
SpectroWe4j list, as compiled and improved by A. Lobel 
(iLobell [20081 . l20l~il ). Where the oscillator strengths (log gf) 
for the same line differed between the lists by more than 
0.25 dex, the value was used that provided an abundance 
best matching the rest of the lines of the same ion. 

To study the Fel/Fen balance, the Fe abundance was 
computed with MOOGIO for three best values of T e g (5750, 
6000, and 6250 K), a range of logg =0.5-3.0 with a step of 
0.5 dex, and a range of micro-turbulences Ltur=3-8 km s -1 
with a step of 1 km s -1 . Only Fel lines with EW^llO mA 
and Fell lines with EW^200 mA were retained for further 
analysis, as they lie on the linear part of the curve of growth. 


3 The following units are used throughout the paper: K for the ef¬ 
fective temperature T e ff, dex[cm s —2 ] for the logarithm of surface 
gravity logg, km s —1 for the micro-turbulent velocity Vtur, and 
dex for the logarithm of the elemental abundance with respect to 
the hydrogen abundance (using designation (X/H) for log e(X) 
on the scale where log e(H) = 12, and [X/H] when expressing 
difference with the Sun: [X/H]=(X/H)—(X/H)©). 

4 http://spectra.freeshell.org/spectroweb.html 
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eW, mA 

Figure 1. Ionization balance for Fe lines, that was used to de¬ 
rive the following atmospheric parameters for IRAS 19135+3937: 
Vtur=5 km s- 1 , log g=1.0, log e(Fe) = 6.50. The plot is shown for 
T eff =6000 K, which is the best-fit value from the fit to the Hy¬ 
drogen lines. Straight lines designate linear fits to the data. 

For each model, the derived abundances from the individual 
lines have been plotted against the EWs. We obtained differ¬ 
ent values for the micro-turbulence when trying to remove 
the abundance trends with E Ws for the neutr al and ionic 
transitions. As discussed in iGorlova et al.1 (l2012h . at these 
temperatures and gravities strong Fel lines are susceptible 
to deviations from local thermodynamic equilibrium (LTE), 
therefore, Vtur from the Fell lines was adopted. The sur¬ 
face gravity is then established by requiring that the average 
Fe abundance from the Fe II lines matches the extrapolated 
to EW=0 abundance from the Fel lines (Fig. [!}. Follow¬ 
ing this procedure for each of the three values of T e s ob¬ 
tained from the hydrogen line analysis, we obtained the fol¬ 
lowing Teff/log <?/E t ur/[Fe/H] solutions: 6000/1.0/5.0/—0.97 
(best model), 6250/1.5/6.0/-0.84, and 5750/0.5/5.0/-1.12, 
where the Fe abundance is expressed relative to the solar 
value of 7.47 in the scale log e(H) = 12. 

Similar rules were used for the rest of the elements to 
combine abundances from the individual lines: for ions the 
abundances of all lines with EW < 200 mA were averaged; 
for neutrals an extrapolation to EW=0 was performed for 
lines with EW <110 mA when possible, otherwise, abun¬ 
dances of lines with EW < 50 mA were averaged, or, if 
such lines did not exist (Col, Znl), lines with EW up to 
110 mA have been averaged. These final abundances for the 
best atmospheric model are given in Table [T] along with the 
difference with abundances obtained using two other closest 
models. 

All elements in IRAS 19135+3937 show sub-solar abun¬ 
dance, ranging from —0.5 for the CNO group to —1.5 for 
some heavy elements. The under-abundance appears to cor¬ 
relate with the condensation temperature of the element, as 
shown in Fig. [2] This is a common phenomenon in pAGB 
disc systems. It is explained by contamination of the photo¬ 
sphere by the re-processed gas from the disc. The gas is poor 
in refractory elements because they condensed into grains. 



0 500 1000 1500 


Figure 2. Abundances of IRAS 19135+3937 (computed with 
the best-fit model T e ff/log3/Vt ur =6000/l.0/5.0) versus conden¬ 
sation temperature. The errorbars represent the line-to-line scat¬ 
ter, while the circle size designates the reliability flag: the larger 
the circle, the more reliable the flag. 

iRao fe Giridhail ihoidh carried an independent analysis of 
the IRAS 19135+3937 composition, in the framework of 
their survey of objects occupying the RV Tau box on the 
IRAS colour-colour diagram. The spectrum that they inves¬ 
tigated was obtained in the phase of line doubling, and so is 
less favourable for measuring EWs than ours. Nevertheless, 
they obtained very similar atmospheric parameters includ¬ 
ing metallicity: T e ff/logg/Kur/[Fe/H]=6000/0.5/4.1/—1.04. 
The depletion pattern that we fi nd for this star is rela tively 
weak, so it is not surprising that lRao fc Giridhail (l2014h did 
not manage to detect it. They did point out the deficiency 
of the a-elements Ca and Ti, that have a higher condensa¬ 
tion temperature than Fe, but were reluctant to ascribe it to 
the depletion effect due to a lack of a stronger Sc deficiency 
and of the Zn enrichment. Our values for these elements are 
more consistent with the depletion pattern, except for Zn, 
but the latter is only represented by 1 or 2 lines in both 
studies. 

The range of temperatures and gravities obtained for 
IRAS 19135+3937 near maximum ligh t allows us to esti¬ 
mate its spectral type fSpTh [Kovtvukh (2007) adapted the 
traditional method of the spectral type determination, that 
rests on the relative strengths of certain lines, to the echelle 
spectra of FGK supergiants. The S/N of our spectrum is not 
sufficient to apply this method directly, but we can use a 
tabulated SpT-T e ft relationship from that study to estimate 
the range of SpTs for IRAS 19135+3937 based on its T e g 
: F6 (6268 K) - F9 (5752 K). In addition, we searched the 
UVES Paranal Observatory Project (UVES POP) archi vcQ 
(Bag nulo et al.l l2003l l for bright field stars w i th sp ectral 
classification adopted from lBuscombe fc Foster! (119951 1. that 

5 http://www.eso.org/sci/observing/tools/uvespop.html 
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Table 1. Chemical composition of IRAS 19135+3937 


z 

Ion 

log e 

[X/H] 

rms 

A [X/H] 

A B 

N 

Flag 

6 

Ci 

8.10 

-0.42 

0.11 

-0.03 

+0.03 

8 

2 

7 

Ni 

7.76 

-0.17 

0.11 

+0.03 

-0.05 

4 

2 

8 

Oi 

8.21 

-0.68 

0.21 

-0.09 

+0.09 

2 

3 

11 

Nai 

5.72 

-0.60 

- 

-0.08 

+0.08 

1 

3 

12 

Mg i 

6.88 

-0.70 

- 

-0.05 

+0.03 

1 

3 

13 

All 

5.60 

-0.87 

0.10 

-0.05 

+0.04 

2 

3 

14 

Sii 

7.06 

-0.48 

0.12 

-0.07 

+0.07 

13 

2 

14 

Sin 

7.26 

-0.36 

- 

-0.01 

-0.09 

1 

3 

16 

Si 

6.57 

-0.55 

0.11 

-0.09 

+0.09 

2 

3 

20 

Cai 

5.20 

-1.16 

0.17 

-0.11 

+0.09 

12 

2 

21 

Sc n 

1.93 

-1.22 

0.14 

-0.26 

+0.22 

11 

1 

22 

Tii 

3.77 

-1.24 

0.01 

-0.21 

+0.19 

2 

3 

22 

Tin 

3.81 

-1.20 

0.12 

-0.24 

+0.19 

24 

1 

23 

Vn 

3.17 

-0.80 

- 

-0.22 

+0.19 

1 

3 

24 

Cri 

4.68 

-0.93 

0.12 

-0.11 

+0.11 

8 

2 

24 

Cr ii 

4.63 

-0.96 

0.07 

-0.16 

+0.13 

13 

1 

25 

Mm 

4.38 

-1.0 

- 

-0.14 

+0.13 

1 

3 

26 

Fei 

6.46 

-1.01 

0.20 

-0.13 

+0.13 

142 

2 

26 

Fen 

6.54 

-0.94 

0.12 

-0.18 

+0.14 

43 

1 

27 

Co i 

3.77 

-1.12 

- 

-0.12 

+0.12 

1 

3 

28 

Nil 

5.35 

-0.88 

0.14 

-0.13 

+0.14 

27 

2 

29 

Cui 

2.83 

-1.38 

- 

-0.2 

+0.2 

1 

3 

30 

Zm 

3.55 

-1.02 

- 

-0.15 

+0.15 

1 

3 

39 

Y ii 

0.66 

-1.56 

0.11 

-0.26 

+0.25 

6 

1 

40 

Zr ii 

1.45 

-1.13 

- 

-0.24 

+0.25 

1 

3 

56 

Ban 

0.92 

-1.21 

- 

-0.3 

+0.26 

1 

3 

57 

La ii 

0.17 

-1.02 

- 

-0.3 

+0.3 

1 

3 

58 

Ceil 

0.33 

-1.20 

0.04 

-0.29 

+0.29 

3 

1 

60 

Nd ii 

0.46 

-1.10 

0.11 

-0.31 

+0.32 

4 

1 

62 

Smu 

-0.04 

-1.03 

0.07 

-0.31 

+0.32 

2 

3 

63 

Eun 

-0.49 

-0.99 

0.06 

-0.26 

+0.27 

2 

3 


Column [X/H] gives abundances for the best-fit model 
-^eff /1°S p'/Vtur— 6000/1.0/5.0, while the A [X/H] column shows 
the response of the abundances to the change in the adopted T e ff, 
with cases A and B corresponding to the models 5750/0.5/5.0 and 
6250/1.5/6.0. The flags in the last column have been introduced 
as follows: 1 - most reliable abundances, for ions where three and 
more lines with EW+200 mA were available for averaging; 2 - 
less reliable abundances, obtained for neutrals by extrapolating 
to EW=0 the abundances from lines with EW^llO mA, or if 
not possible, by averaging lines with EW^50 mA; 3 - least reli¬ 
able abundances, obtained by averaging less than three lines with 
EW^200 mA and EW^llO mA for ions and neutrals, respec¬ 
tively. The number of used lines is indicated in the last but one 
column. The rms column marks the mean deviation either from 
the average value or from the interpolated line. 


would be similar to IRAS 19135+3937. Given the peculiar 
abundance pattern of IRAS 19135+3937 and a wide range 
of line widths exhibited by supergiants, we did not intend 
to find a precisely matching standard, but rather wanted to 
independently verify SpT of IRAS 19135+3937 implied by 
its T e fj. The spectrum of HD 108968 (F7Ib/II) appeared to 
be the closest match to our spectrum of IRAS 19135+3937, 
hence we adopted SpT F7±2 I/II for the latter. 


3.2 IR excess 

The SED of IRAS 19135+3937 is shown in Fig. [3] Besides 
our measurements shown for the maximum and minimum 
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Figure 3. Spectral energy distribution based on our photom¬ 
etry and the measurements from the VizieR database. Dashed 
line: Kurucz photospheric model with our spectroscopically de¬ 
termined parameters T e fj/log g/[M/H]=6000/1.0/—1.0, reddened 
by FA B — V) = 0.18, and shifted to match V max of our epoch of 
observations. 


light, we also plot photometry from the VizieR database. 
In particular, the measurements in the Johnson’s B and V 
pass-bands originate from the "A ll -sky com piled catalogue 
of 2.5 million stars” (iKharchenko fe Roeserl 120090 . Another 
measurement in the Johnson’s V, and a measurement in 
the Cousin’s I pass-band originate from the “The Amateur 
Sky Surve y (TASS ) Mark IV patches photometric catalog, 
version 2” jDroege et al .1 [jQOfii) . The rest of the data points 
arise from the annotated ground and space-based missions. 
The scatter in the optical is due to the variability of the 
source. 

Using the photospheric parameters determined in Sec. 
13+1) , we can estimate the stellar contribution in the SED and 
identify any flux excess. The photosphe ric contr i bution in 
Fig. [3] is represented by a Kurucz model llCastelli fc Kuruc3 
120031 with T e ff = 6000 K, logp = 1.0, [M/H]= —1.0, and 
Vtur = 2.0. The model was reddened by E(B — V) = 0.18, 
which was deduced from the comparison of the observed 
B — V colour in the maximum light with the expected colour 
for an F7 sup ergian t (Sec. []}. We used reddening law of 
ICardelli et ajj j[l989[l with Ay/E{B — V) = 3.1 with modi¬ 
fication of O’Donnell 11994) in the optical to NIR regime. It 
is clear that IRAS 19135+3937 has a strong IR excess. The 
SED is typical for pAGB stars with dusty discs, where the 
excess usually starts at 2 pm, peaks at around 10 pm and 
drops with a black-body slop e towards the far-IR, indicating 
the presence of large grains llDe Ruvter et al1l2005h . 


4 PHOTOMETRIC VARIABILITY 

In Fig. U we show our light and colour curves of 
IRAS 19135+3937, that are clearly variable. The variations 
in all three bands are periodic and occur in phase. Using 
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JD-2400000 

Figure 4. U-band light, and U — B, B — V colour curves for 
IRAS 19135+3937 using our photometry from 2012—2013. 
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Figure 5 . U — B and B — V colour variations with the U-band 
magnitude for IRAS 19135+3937 in 2012-2013. 

the PeriodO^I program of iLenz fe Bregeil 120051 ). that re¬ 
lies on the discrete Fourier transform for period determina¬ 
tion, we obtained P = 126.53 ± 0.40 d for all three pass- 
bands. The peak-to-peak amplitude increases with wave¬ 
length: AU = 1.0 mag, A B = 1.1 mag, AV = 1.2 mag. 
The average magnitudes in 2012-2013 were: < U >= 12.17 
mag, < B >= 11.80 mag, < V >= 11.23 mag. The rela¬ 
tionship between brightness and colour is striking (Fig. [5]): 
when the star fades, it gets bluer. Clearly, this behaviour can 
not be explained simply by variations in the extinction, as 
it would produce the opposite effect. 

To verify how IRAS 19135+3937 behaved in the past, 
we searched for time series photometry in the publicly avail¬ 
able catalogues. In 1989-1992 it w as observed by Hippar- 
cos. In Tycho-2 llHpg fc et al.ll2000l ) it has the designation 
TYC 3125-2395-1. One hundred and eighty measurements 
are given, as well as the mean magnitudes < Bt >= 10.886 
and < Vt >= 10.200. Convertin g to th e the Johnson sys¬ 
tem using the transformation of iBesse 2 1 I 2000 I) , we obtain: 
< B >= 10.75 mag and < V >= 10.13 mag. The errors of 

6 http://www.univie.ac.at/tops/Period04/ 


the Tycho measurements, however, are too large (av T — 0-29 
mag and ob t = 0.32 mag) for a reliable period determina¬ 
tion. 

In 2003-2004 IRAS 19135+3937 was monitored by 
TASS. Forty-six measurements in Johnson-Cousins V and 
Ic p ass-bands for the first tim e revealed variability of the 
star llSallman fc Droegel l2002 ). We obtained the following 
values of the period for this data set: 127.11±0.95 d (17- 
band) and 127.11±0.71 d (/c-band). According to these ob¬ 
servations, the V — Ic colour did not exhibit sinusoidal vari¬ 
ations with amplitudes larger than the uncertainties of the 
photometric observations (±0.07 mag), in contrast to the 
pronounced U — B and B — V colour variations observed by 
us in 2012-2013. 

Between May 2006 and January 2008 
IRAS 19135+3937 was observed by ASAS. The mea¬ 
surements are given in “The catalogue of variable stars in 
the Kepler field of view’fj]. The source with a designation 
191512+3942.8 has a period P=128.8 d, mean magnitudes 
< V >= 11.269 mag, < I >= 10.441 mag, and the vari¬ 
ability type QPER. (a semi-regular variable with a stable 
period). The observations in both bands unfortunately 
are not simultaneous, and the behaviour of the V — I 
colour with phase depends on the details of interpolation. 
Nevertheless, the V — I variations did not exceed ±0.04 
mag, in agreement with the TASS observations. Using the 
latest data for the smallest aperture size of 30”, which is 
similar to ours and avoids contamination from a nearby 
star, we obtained the period values of 122.06±0.65 d for 
the U-band and 127.15±1.06 d for the /-band. One should 
note that the quality of the U-band measurements in the 
ASAS catalogue is very poor - 91% of data points have 
the lowest grade ’D’ and only 8% grade ’A’. The /-band 
measurements are much better, with 67% having grade ’A’. 
Hence, we adopt P = 127.15 d for this data set. 

In Fig. [6] we show 4 years of Kepler photometry of 
IRAS 19135+3937 (= Kepler 4644922) carried out by the 
Kepler satellite. We plot 18 long-cadence data sets that were 
available in the STScI archive at the beginning of 2014. Due 
to the fact that the time-scale of variations exceeds the du¬ 
ration of a Kepler observing quarter, the de-trending proce¬ 
dure performed by the Kepler pipeline is invalid. Raw fluxes 
are therefore plotted, where small offsets between adjacent 
data sets result from the imperfect calibration. Applying 
Period04 to these data, we obtain P = 127.497 ± 0.014 d. 
The Kepler light curve is the most precise of all observations. 
It is extremely smooth and shows that the brightness varia¬ 
tions are very regular, but the shape is not fully repeatable 
from cycle to cycle. 

The summary of all photometric observations (except 
Tycho-2) is given in Table [5] As can be seen, there is no 
indication from this data that the period changed over the 
past 10 years. The shape of the light-curve and the mean 
magnitude, however, did change from cycle to cycle, con¬ 
firming the semi-regular classification of the star. In Fig. [7] 
we plot all observations in the U-band together. Tycho-2 
measurements have been omitted from the plot due to the 
large error-bars and the Kepler ones due to a non-standard 
pass-band. We performed a period search on this combined 


' http://www.astro.uni.wroc.pl/ldb/asas/kepler.html 
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Figure 6. Kepler light curve. Colours mark separate acquisition 
sequences. 



Figure 7. Combined V-band light curve based on the TASS, 
AS AS and our photoelectric photometry (SAI). 

data set (200 data points), by preliminary scaling each of the 
three data sets to their average magnitude, and obtained 
P = 127.04 ± 2.6 d, which, as will be shown in the next 
section, coincides with the spectroscopic period. 

Our simultaneous UBV photometry enabled us to 
study colour variations in IRAS 19135+3937 for the first 
time. Fig. [8] shows an area on the U — B, B — V two-colour 
diagram traversed by IRAS 19135+3937. T he sequences of 
superg iants and bright giants according to ISchmidt-Kaleii 
d 1982li are drawn for reference. From our spectroscopic anal- 

Table 2. Summary of the photometric observations 



B-V 


Figure 8. Location of IRAS 19135+3937 on the U — B vs. B — 
V diagram, de-reddened by E(B V) — 0.18 and fitted with 
a straight line to guide the eye. The corresponding reddening 
vector is shown as a red line, following equation E(U — B)/E(B — 
V) = 0.72 + 0.05 E(B — V). Black solid and dashed lines designate 
sequences of supergiants (lb) and bright giants (II), respectively. 

ysis it follows that at maximum light the star should be of 
spectral type F7 (±2). Accord ing to the superg i ant ca libra- 
tions of ISchmidt-Kalerl d 198 2t i and iKovtvukhl l|2007l ). this 
corresponds to ( B — V)o = 0.48 ±0.1 mag. The observed 
(B — V) at maximum light is ~0.66, hence E(B — V) = 
0.18 ± 0.1 mag. Assuming that this reddening is predom¬ 
inantly interstellar 2 _andado£ting the standard reddening 
law of iHiltner fe .Tohnsonl (1 19561) : E(U - B)/E(B - V) = 
0.72 + 0.05 E(B — V),- we have de-reddened all points in 
Fig. [8] by this value. As can be seen, the colours vary consid¬ 
erably throughout the cycle, but follow neither reddening 
nor temperature sequences. While in the maximum light 
both colours can be explained by a moderately reddened 
F7 lb supergiant, in the minimum light the star can not be 
brought to the same locus on the UBV diagram even when 
allowing for a different amount of reddening. Pulsations can 
hardly explain this behaviour either, as such stars normally 
get redder at minimum light. Could the blue colours at min¬ 
imum light be due to the contribution of an early-type com¬ 
panion? In the next section we will examine our time series 
of spectra to investigate this possibility. 


survey 

P 

< V > 

AV 

Tstari^ 

AT 


(d) 

(mag) 

(mag) 

(d) 

(d) 

TASS 

127.11+0.95 

10.69 

0.9 

52782 

477 

ASAS 

127.15+1.06 

11.36 

0.7 

53884 

598 

Kepler 

127.50+0.01 

N/A 

N/A 

54953 

1471 

SAI 

126.53+0.40 

11.23 

1.2 

56093 

519 


a First day of observation, JD-2400000 


5 SPECTROSCOPIC VARIABILITY 
5.1 Radial velocities 

To determine RV of IRAS 19135+3937, we used a cross¬ 
correlation method with a mask containing lines of a G2 star. 
The correlation was performed on the order-by-order basis, 
after which the cross-correlation functions (CCFs) from the 
central 20 orders were merged into one to represent an av- 
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Figure 9. Cross-correlation functions arranged according to the 
orbital phase, with time running downwards. Throughout the 
paper the phases have been defined as follows: 0 = (JD — 
2454980.42464)/127.08, with 0 = 0 corresponding to RVmax . The 
colours on the right panel designate the year of observation to il¬ 
lustrate the cycle-to-cycle variations. The solid and dashed lines 
trace the main and the secondary CCF components, respectively. 


erage line profile. Normally, this profile would be fit with a 
Gaussian to obtain an RV measurement. However, in many 
phases the CCFs turned out to be of a more complex shape. 
While some CCFs have the form of a single absorption com¬ 
ponent, either strong and narrow or shallow and broad, in 
other CCFs two separate components can be clearly seen. 
The determination of the CCF centroid is therefore not 
straightforward. 

To investigate whether this behaviour could be peri¬ 
odic, we first measured a flux-averaged value of the RV for 
each CCF, using (1 — F\) 2 as a weighting function. Apply¬ 
ing Period04 to these data revealed several peaks on the 
periodogram, of which the first two strongest ones are at 63 
d and 129 days. The latter coincides with the photometric 
period, while the former reflects the fact that the lines split 
twice per photometric period. Arranging the CCFs according 
to the phase calculated with the longer period allowed us to 
identify the main component (which is usually the stronger 
one) and the secondary component in the split CCFs (Fig. 
®. We then re-determined RVs separately for each compo¬ 
nent using a double-Gaussian fit. In several cases it was not 
possible to identify the secondary component, because it was 
in the noise or merged with the main component. In such 
cases a single-Gaussian fit was performed, and the RV was 
assigned to the main component. The resulting RVs are given 
in Table A2. 

In order to determine whether the main RV component 
could be due to the orbital motion of the brighter star (the 
primary) in a binary system, we fit it with a Keplerian orbit, 
which included a re-evaluation of the period. The fit is very 


Table 3. Orbital elements for the primary (pAGB) component 
of IRAS 19135+3937. 


Parameter 

Value 

(7 

p (d) 

127.08 

0.08 

a sini (AU) 

0.20 

0.007 

/M (M 0 ) 

0.07 

0.007 

K (km s -1 ) 

17.71 

0.5 

e 

0.14 

0.02 

w (°) 

66 

14 

T 0 (JD) 

2 454 997.8 

6.8 

7 (km s -1 ) 

1.3 

0.4 

^red 

1.11 


R 2 

95.4% 



Listed are orbital parameters with their uncertainties, reduced 
chi-square, and the coefficient of determination. 

good - see Tableland Fig. [TO] The larger residuals from the 
fit between phases 0.4 - 0.9 are due to the broader CCFs. 
To estimate the uncertainties on the orbital parameters, we 
carried out 1000 Monte-Carlo realisations to simulate the 
observed RVs. For each date we used a Gaussian distribution 
centred on the observed RV and a adopted as follows: 0.7 
km s _1 for the orbital phases between 0 = 0.95 — 1.35 and 
3.5 km s~ x for 0 = 0.35 — 0.95, based on the residuals from 
the Keplerian fit. For each realization the orbital parameters 
were determined, and the a of the distribution of a given 
parameter was adopted as the parameter’s uncertainty. Fig. 
d illustrates this procedure for the period and eccentricity. 

The RV curve of the secondary CCF component is much 
less certain, in particularly between phases 0 — 0.3, so we did 
not attempt to model it. However, we can use the deduced 
mass function from the fit to the RV of the main component 
to estimate the mass of the putative companion. Adopting 
0.5 Mg for the mass of the pAGB primary (a mass of a 
typical white dwarf), with the mass function of 0.07 Mg 
we obtain masses of the companion between 1.1 and 0.4 
Mg, depending on the adopted inclination of 30° and > 
65°, respectively. The companion has a comparable or larger 
mass than the primary. 

IRAS 19135+3937 turned out to be yet another disc 
system that is a binary. The line behaviour is reminiscent 
of the double-lined spectroscopic binaries, but not identical: 
the relative strengths of the two CCF components in our 
case vary, both with phase and from cycle to cycle. While the 
main component can be safely attributed to the pAGB star, 
the association of the secondary component with a physical 
companion is not straightforward. 

In Fig. [T^] we overlay RVs of the main CCF component 
and the Kepler photometry, that were obtained over the 
same time-span. Despite the cycle-to-cycle variations in the 
light curve, an offset of a quarter of a period can be clearly 
seen with the RV curve: the superior conjunction of the pri¬ 
mary star (0 ~ 0.2) coincides with the maximum light, while 
the inferior conjunction (0 ~ 0.7) with the minimum light. 
Brightness declines are therefore consistent with the obscu¬ 
ration of the primary by the inner disc wall. 

Over half of the orbit centred on the minimum light 
(0 = 0.4 — 1.0) one also observes an increase in the RV 
scatter. As can be seen in Fig. 0 this is not due to the 
decreased S/N, but due to the fact that the secondary com- 
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Figure 10. Radial velocities for the main component of the CCF 
(orange circles) and for the secondary component (blue). The for¬ 
mer were fit with a Keplerian orbit, which is drawn in the upper 
panel, while the residuals from it are plotted in the lower panel. 
The size of the circles is proportional to the component’s strength. 
Horizontal lines mark d=3cr deviation from the fit orbit. 



Period, P + 127 (d) 


Figure 11. Distribution of values for eccentricity and orbital pe¬ 
riod based on 1000 Monte-Carlo simulations of the RVs of the 
main CCF component. Solid lines mark a Keplerian solution ob¬ 
tained from the observed RVs, whereas the dotted lines denote a 
simulated RV set that resulted in the smallest \ 2 value. 


ponent of the CCF becomes of comparable strength to the 
main component, resulting in the shallower, smeared com¬ 
bined profiles that are difficult to fit. This is not surprising 
considering that at </> = 0.7 the primary becomes obscured, 
and the spectrum of the companion should become more vis¬ 
ible. In Sec. E3 however, we will present some arguments 
that hinder a definite identification of the secondary CCF 
component with the spectrum of the companion. 



Figure 12. Radial velocities (orange circles) and Kepler fluxes 
(red to blue points) plotted against the orbital phase. The colour 
scheme is the same as in Figs.[9land ll4l The horizontal line marks 
the systemic velocity. 


5.2 Emission lines 

There are a number of variable emission lines in the spec¬ 
trum of IRAS 19135+3937. Fig. [13] presents an overview 
of the discussed features in the form of dynamic spectra. 
The upper panel shows the behaviour of the photospheric 
lines using the example of the cross-correlation function 
and one strong Ban line; the middle panel shows emission- 
absorption profiles of Ha and one component of the NaD 
doublet; the bottom panel compares static emission lines of 
Fe and TiO with the nearby dynamic photospheric lines. 

The most prominent emission line is Ha (higher mem¬ 
bers of the Balmer series show a similar behaviour, but to 
a much smaller degree and confined to the line core). Ha 
profiles, arranged according to the orbital phase, are shown 
in Fig. Q3] Over one half of the period, centred on (j> = 0.7 
(pAGB inferior conjunction), the profiles take a shape of a 
double-peak emission, with a distance between the peaks of 
~100 km s^ 1 and a width at the base of at least 300 km s -1 
(before photospheric absorption subtraction). Alternatively, 
such profiles can be interpreted by a broad emission profile 
with a narrower superimposed absorption component. This 
absorption can not be photospheric due to the fact that it 
does not follow the orbital motion, but is permanently 13 
km s~ x blue-shifted relative to the systemic velocity. Over 
another half of the period the emission in Ha is overtaken by 
a broader, asymmetric absorption, with a blue wing that ex¬ 
tends to over 350 km s” 1 near <f> = 0.2 (pAGB superior con¬ 
junction). This b eha viour is identical to that in BD+46°442 
llGorlova et al.ll2012l ). 

In some spectra we also spotted weak emission in a num¬ 
ber of low-excitation metallic and molecular lines. To study a 
possible correlation with the orbital phase, we first averaged 
spectra within each of 10 phase bins to increase the S/N. 
These combined spectra arranged according to the phase 
are shown in Fig. 1151 Three examples of emission features 
are given: the resonance lines of Nal, low-excitation Fel line 
8047.62 A (xiow — X-up — 0.9 — 2.4 eV), and three close TiO 
band-heads. 
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Figure 13. Dynamic spectra of the representative features in the spectra of IRAS 19135+3937: top panel shows the behaviour of the 
photospheric absorption lines, middle panel - lines dominated by the circumstellar emission and absorption, bottom panel - pure emission 
lines (with the neighbouring photospheric lines of Fel and Cl shown for reference). Colours represent continuum-normalized fluxes, with 
unity corresponding to the continuum level. Spectra are arranged according to the orbital phase with period 127.08 d; one period is 
shown twice; time runs down. The following lines are drawn to guide the eye: vertical line for the centre-of-mass velocity, solid curve - 
for the Kepler’s solution of the main CCF component, dashed curve - of the secondary component. Short horizontal dashes designate 
observed phases. 
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Figure 14. Ha as a function of the orbital phase. Vertical line 
marks the systemic velocity. 


The Na D lines in IRAS 19135+3937 are dominated by 
a composite “shell” absorption component, which is common 
for pAGB stars. The photospheric component is only visible 
around <j> = 0, when it is maximally redshifted. Through¬ 
out the rest of the orbit it is hidden in a cluster of narrow 
absorption lines, that are permanently blue-shifted relative 
to the systemic velocity. Different features can be better 
traced in the dynamic spectra, as shown in Fig. [13] (middle 
right panel). There, the Na D photospheric absorption can 
be traced up to <j> = 0.4, which corresponds to the least ob¬ 
scured part of the primary’s orbit, when other photospheric 
lines are strongest as well. At least two non-photospheric ab¬ 
sorption components can be disentangled: at —38 and —20 
km s _1 , that can be both of inter- and circumstellar origin. 


What is not so common is the appearance of the broad emis¬ 
sion wings between <j> = 0.6 —0.8 (in Fig |13l thev can be seen 
as a flux excess over the mean profile), that are reminiscent 
of emission in Ha. 

Some weak emission lines are found in the red part of 
the spectrum, where they better stand out against the drop¬ 
ping continuum. The three most notable lines can be asso¬ 
ciated with two close low-excitation multiplets (x«p = 2 — 3 
eV) of Fel, with rest wavelengths of 6400.32, 6498.94, and 
8047.62 A. Many more can be hidden in the telluric fea¬ 
tures and noise. Careful examination reveals that they have 
double-peak, but relatively narrow, profiles (the distance be¬ 
tween the peaks is ~ 30 km s -1 ), that are permanently cen¬ 
tred on the systemic velocity (within 3 km s _1 ). The latter 
two lines were also noted by us in BD+46°442, but in the 
case of IRAS 19135+3937 they are stronger and clearly vari¬ 
able. In the continuum normalized spectra, they appear to 
vary in strength with the photometric/orbital period. They 
are 2.5 times stronger in the minimum light than in the max¬ 
imum, which coincides with a one magnitude amplitude of 
the photometric variations. This means that the line emis¬ 
sion flux is constant, and the apparent line variability is due 
to the variations in the continuum flux. 

Finally, we detected weak emission from three TiO 
band-heads belonging to the 7(0,0) system. Similarly to 
the atomic emission lines, they become strongest during the 
minimum light, and appear to be stable in the RV. The in¬ 
dividual transitions, however, are blended with each other, 
so it is not possible to gauge the kinematics of the emitting 
region. 


6 DISCUSSION 

We showed that IRAS 19135+3937 possesses typical charac¬ 
teristics of a pAGB star with a dusty disc: it is a low-gravity 
F star with large IR excess but low reddening, it is deficient 
in refractory elements, and it has a strong and variable Ha 
profile indicative of mass loss. In particular, it is ver y similar 
to BD+46°442 described in our earlier work (IGorlova et ali 
1201^1 . The only substantial difference between the two is 
that IRAS 19135+3937 is a recognized variable star due 
to the larger amplitude of the light curve. The HERMES 
survey uncovered periodic, but complex spectral behaviour 
in both stars. In the following subsections we will elabo¬ 
rate upon our model of an interacting binary system that 
we developed for BD+46°442 to explain additional aspects 
exibited by IRAS 19135+3937: brightness variations (Sec. 
nm . stationary emission lines (Sec. 16.21) . photospheric line 
splitting and peculiar colour variations (Sec. 16.31) . 

6.1 Binarity and SRd variability 

For the first time we have demonstrated that 
IRAS 19135+3937 is also variable in radial velocity, 
with a period coinciding with the photometric period. We 
can now use the same argument as for BD+46°442 to 
disregard pulsations as the cause of variability. Integrating 
the RV curve of the main CCF component over half of 
the period, we obtain a displacement ~90 Rq, which is 
comparable with the radius of the pAGB star itself and 
would be prohibitively large for pulsations. And if line 
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Figure 15. Emission lines as a function of orbital phase. Black: an average of all spectra falling within a given phase bin; red: an average 
spectrum over the entire orbit. All spectra have been normalized to unity at the continuum and shifted vertically according to the orbital 
phase. Vertical lines mark wavelengths for the identified features at the systemic velocity. 


splitting were due to the propagation of shock waves, as in 
RV Tau stars, one would have to explain why this star does 
not show other characteristic features of these variables, 
such as interchanging deep and shallow minima, and, most 
notably, variations in the effective temperature. 


The constancy of T e g is demonstrated in Fig. 1161 where 
we compare a part of IRAS 19135+3937 spectrum in the 
phases of minimum and maximum light with the spectra of 
two standard F supergiants fro m the UVES POP archive: 
HD 74180 (F3 la , Malarodalli 975)1 ) and HD 108968 (F7 lb/II, 
iHouk fc Cowlevlll975l f. The standards were chosen to have 
SpTs expected for IRAS 19135+3937 based on its B — V 
colour in those phases. It can be seen that in the late-F stan¬ 
dard the lines of neutral species become noticeably stronger 
than in the early-F one. In the spectra of IRAS 19135+3937, 
however, this effect is not observed, hence, the tempera¬ 
ture must have remained constant between the two op- 
posite phases. Another confirmation of this fact is that 
iRao Sz Giridhail (l2014tl studied IRAS 19135+3937 in the 


intermediate phase (</> = 0.49) and obtained precisely the 
same value of T e g as us. 


We conclude that the RV variations in 
IRAS 19135+3937 are due to the orbital motion of 
the giant primary around a much fainter companion, 
which is very common among pAGB stars with discs. 
The periodic semi-regular dimmings can be explained by 
the obscuration of the primary by the circumbinary disc 
matter, most likely by the puffed-up inner rim during 



A, A 


Figure 16. Comparison of IRAS 19135+3937 spectra in the 
phases of minimum and maximum light with each other and with 
two spectral standards. 
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the inferior conjunction. The absence of a flat part on 
the light curve near maximum light could be due to the 
permanent obscuration or presence of a large amount of 
scattered light. This explanation has been already proposed 
for a few other long-period pAGB variables with the same 
phase shift between the lig ht and RV curves, such as HR 
4049 (IWaelkens et all ~1991f). HD 52961 jyan_Winckel et all 


1 199 fill . and EN TrA (I Van Winckel &: et alj 2009h . The very 
smooth Kepler light curve indicates that the photosphere 
of IRAS 19135+3937 is very stable. 


6.2 Accretion onto companion and gas discs 

Furthermore, a number of spectroscopic features in 
IRAS 19135+3937 point to the presence of an active mass 
transfer between the primary and companion. In particular, 
as was described in Sec. E2 near <j> = 0.25 Ha develops 
a spectacular P Cyg-like profile. With our phase conven¬ 
tion, where phases are counted from the time of the maxi¬ 
mum RV of the primary, <j> = 0.25 for a nearly-circular orbit 
corresponds to the giant’s superior conjunction (where it is 
farthest aw ay from us, whi l e the puta tive companion is in - 
between). In lGorlova et al.l d2012l l and lGorlova et ail] J 20131 ) 
we describe a few other systems from the HERMES survey 
and from the l iterature w ith a sim ilar behaviour of Ha. As 
was shown bv lThomas et all d2013h for the central star of the 
RR nebula, this phenomenon can be explained by a wide- 
angle jet originating at the secondary. The lobe pointing in 
our direction will be periodically projected against the gi¬ 
ant primary and produce the observed transient blue-shifted 
absorption. The jet is likely powered by accretion from the 
pAGB primary to the secondary. 

The interpretation of the double-peaked emission lines 
is not so straightforward. Usin g STIS sp ectro graph on the 
Hubble Space Telescope , iThomas et al.l ll201ll) spatially re¬ 
solved narrow emission (with the distance between the peaks 
of 12 km s -1 ) in the NaD lines of the RR, and deduced that it 
is produced in the distant parts of a bipolar outflow, that are 
seen in the direct light. On t he other hand, the much broader 
emission in Ha according to 1 Witt et ali (2009) could form in 
the parts of the lobes that are closer to the binary. The RR, 
however, is not a typical disc object, and the appearance of 
some features may be affected by its nearly perfect edge-on 
orientation. 

In Figure [17] we show another possibility to explain the 
formation of th e double-peaked emission profiles - in a Kep- 
lerian disc. ISma kl l|l969h presented a simple kinematic model 
of purely gaseous, optically thin disc of constant thickness 
with density linearly dropping to zero at the outer edge. 
Using this formulation, we could successfully fit Fe emission 
line profiles in IRAS 19135+3937 with the following disc pa¬ 
rameters: the ratio of the inner to outer radius of 0.15, and 
the velocity at the outer edge of 12 km s” 1 . Depending on 
the total mass of the system (0.9-1.6 M@) and the inclina¬ 
tion angle (30-70°), this results in a disc with Ri n = 0.2—1.3 
AU and R ou t = 1.4 — 8.7 AU. 

Is this Fe disc circumcompanion or circumbinary? In 
Sec. 15.11 we deduced that the semi-major axis of the pAGB 
primary is 0.2 - 0.4 AU and the secondary is likely more 
massive than the primary. Hence, the distance between the 
companions does not exceed 0.4 - 0.8 AU, which is smaller 
than the outer radius of the Fe disc. This fact, together 


with the constancy of the RV of the emission lines, imply 
that this disc is circumbinary. Furthermore , when app lying 
a 2D radiative transfer code (jGielen et al.l i2007l . 2009 1 to 
fit the SED of IRAS 19135+3937, we obtain that the inner 
radius of the circumbinary dusty disc is 5 - 10 AU. Thus, the 
double-peaked Fe emission lines may signal the presence of a 
gaseous Keplerian circumbinary disc, that is nestled within 
the sublimation boundary of the dusty disc. 

In the prototypical disc object 89 Her the above dis¬ 
cussed me tal lines form part of a much richer em ission line 
spectrum dClimenhaga et ah 1987|^ Kippeij jgO 1 il l . Based on 
the constancy of the RV. IWaters et, al.l 1 1993l l were first to 
propose that emission cou l d origi nate in a circumbinary disc. 
Indeed, iBuiarrabal et all (l2007i l possibly detected such a 
disc as an unresolved component in the interferometric maps 
of CO. The existence of gas inside a dusty disc has been long 
anticipated in the framework of the re-accretion hypothesis, 
designed to explain a depletion pattern in some disc hosts. 
The observational evidence of circumbinary gas, however, 
is largely missing. CO studies nor mally probe material o n 
a much larger scale (10 3 -10 5 AU. IBuiarrabal et al.ll2013ll . 
which could have been ejected in the preceding AGB stage. 
In contrast, metal emission lines probe gas on AU scales, 
and therefore provide a better insight into the current mass 
loss/accretion. 

Furthermore, the ~10 times wider emission in Ha may 
indicate the presence of gas well inside the binary’s orbit. 
Giving that there is a jet emanating from the companion, 
this hot gas could form in the circumcompanion accretion 
disc. Interestingly, the rare appearance of TiO in emission 
has been al so associated with t he pr esence of an accretion 
disc. Thus, iHillenbrand et aj i2012i l proposed to explain 
TiO emission in some young stellar objects (YSO) and Be 
stars by evaporating disc material at the base of the outflow, 
where it is lifted up and exposed to the UV radiation from 
the accretion disc. Besides IRAS 19135+3937, TiO emis¬ 
sion has been recently discovered in a few can did ate pAGB 
binaries in the Magellanic Clouds llWood et alJl20l3 ). 

This raises some important questions about the type of 
the accretion (a wind, a Roche-lobe one, or perhaps even ac¬ 
cretion from the circumbinary disc), and whether the three 
discs (the circumcompanion, circumbinary, and the outer 
dust+gas one) could be possibly related to each other. The 
exploration of these possibilities, however, is beyond the 
scope of this paper. 


6.3 Reflection of the circumbinary disc 

What is the nature of the companion and could it explain 
the remaining peculiarities of the system: the secondary set 
of spectral lines and the blueing of the colours in the light 

minima? _ 

We applied the FDBinary cod^3 ( llliiic et alj|2004l ) to 
our spectral data set in the attempt to separate the two 
components. The code performs separation of spectra in 
a spectroscopic double-lined binary star in Fourier space 
without the use of template spectra. The input consists of 
the continuum-normalized spectra, the relative fluxes of the 
two components ("light factors") per observation, that were 


http: //sail.zpf.fer.hr/fdbinary 
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Figure 17. Mean, continuum-subtracted observed profile of an 
Fe emission line in IRAS 19135+3937 (dots) versus the best- 
fit model profile from a Keplerian disc (solid). Both profiles have 
been normalized to unity at maximum flux. The broad absorption 
feature in the observed spectrum is due to the averaging of a 
neighbouring photospheric line over the entire orbit. 

taken from our double-gaussian decomposition of the CCFs, 
and the first approximation for the orbital parameters. In 
the output we obtained two identical sets of lines, whose 
shape and relative depths depended somewhat on the se¬ 
lected regions and the number of iterations, but in no case 
was there an indication that the sets would correspond to 
two different spectral types. This fact, along with the con¬ 
stancy of temperature in the opposite conjunctions, implies 
that the companion must have an identical SpT and com¬ 
parable luminosity to the primary. And yet, only one star 
is seen in the light curve (no substantial secondary mini¬ 
mum/maximum is observed over the RV period). Also from 
the evolutionary point of view it is very unlikely to find a 
system composed of two yellow evolved stars, as pAGB stage 
is extremely short (^1000 years). 

To resolve these difficulties, we propose an alternative 
explanation for the secondary component of spectral lines, in 
IRAS 19135+3937, BD+46°442, RR and potentially other 
inclined disc systems with periodically distorted lines. In 
this picture the companion is undetected, probably a non- 
evolved dim MS star. The secondary set of lines, which is 
most pronounced near the inferior conjunction of the pri¬ 
mary, would be the light of the primary that is reflected 
off the inner wall of the disc and scattered into our line of 
sight. The reflected image will move in the opposite direction 
from the star and disappear during the superior conjunction 
of the latter, explaining the CCF behaviour. 

A very similar phenomenon of line splitting due to 
the presence of reflected light w as described for o ne clas¬ 
sical T Tauri binary KH 15D llHerbst et al.l 120081 ). Sim¬ 
ilar to IRAS 19135+3937, KH 15D experiences periodic 
semi-regular deep declines, that are explained b y the or¬ 
bital motion inside a precessing w arped disc llHerbst et al.l 
l20ld : IWindemuth fe Herbstll2014l ). The SEDs of pre-main- 
sequence (pre-MS) objects are often indistinguishable from 
the pAGB ones, poin ting to the similarity of their discs 
(Ide Ruvter et al]|2006h . One such property is the presence 
of grains that are larger than in the interstellar medium. 


Based on the comprehensive modelling, iHerbst et al.1 d2008l l 
concluded that the reflected light in KH 15D provides the 
best evidence that the particles in the disc grew to at least 
1 mm in size. For pAGB discs the same was previously in¬ 
ferred from the black-body slope of the far-IR excess, which 
can only be measured for brighter sources. Modelling of the 
scattered line profiles may prove to be a useful novel tool 
for studyi ng com p ositio n and kinematics of post-AGB discs 
fe.g. I Grinin et al]l2006h . While in YSOs the particle growth 
is the first step toward the planet formation process, the 
evolution of grains in post-AGB discs is unknown. 

Large particles, however, produce gray scattering and 
settle to the mid-plane of the disc, and hence can not ex¬ 
plain the blueing of the system in the minimum light when 
the primary sinks behind the disc edge. No trace of a hot 
companion is observed in our spectra either. We propose 
that the blueing is due to the scattering by the inner wall of 
the disc. Since we do not resolve the system, the total flux in 
the telescope beam is always a sum of direct and scattered 
light in our line of sight. During the inferior conjunction of 
the primary (phase of minimum light) the amount of bluer, 
back-scattered (by the farthest side of the inner disc wall) 
light may exceed the amount of direct, reddened light from 
the primary, which is possible because, unlike us, the inner 
disk wall always sees the star unobscured. As a result, in 
this phase the system will appear bluer to us than in the 
unobscured phases. 

The colours and the level of the scattered light bring 
stringent constraints on the size distribution and on the 
chemo-physical properties of the dust grains, as well as on 
geometry of the inner dusty disc. The reproduction of these 
observables will need detailed radiative transfer models in 
which the angle and colour-dependent scattering, as well as 
dust settling need to be incorporated. This is outside the 
scope of this paper, however, and will be the subject of 
a subsequent analysis. The high level of optical scattering 
seems to be a common property of the circumbinary discs, 
as also in 89 Her this h a s been detec ted using optical inter¬ 
ferometry. iHillen et al.l (120131 . l2014h spatially resolved the 
optical scattering component and deduced that as much as 
40% of the optical flux of 89 Her is due to scattered light. 

6.4 Circumstellar geometry 

In Fig.[TH]we present an artist’s impression of the circumstel¬ 
lar environment of IRAS 19135+3937 based on the discus¬ 
sion in the previous subsections. Like in BD+46°442, in this 
system the giant primary is transferring mass to a, likely 
unevolved, companion, which results in the production of 
a pair of jets emanating from the secondary. The system 
is surrounded by a dusty disc, and possibly also a smaller 
gaseous disc. The new element in this picture is the high 
level of optical reflection from the primary on the wall of 
the circumbinary disc, that mimics a twin companion in 
the spectra. This reflected light is observed with a different 
Dopplershift and hence can be differentiated from the com¬ 
ponent of direct light. The reflected spectrum and the light 
amplitude are more pronounced in IRAS 19135+3937 than 
in BD+46°442 likely due to the higher inclination of the for¬ 
mer. The contribution of this reflected light is maximal at 
minimum light, which indicates an efficient back-scattering 
and makes the system to appear bluer when fainter. 
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Figure 18. Artist’s impression of IRAS 19135+3937 interacting binary system, depicted in the phase near light minimum. The pAGB 
primary is about to sink behind the dusty disc edge, giving way to the spot of reflected light on the opposite side of the circumbinary 
disc. 


7 CONCLUSIONS 

We presented contemporaneous photometric and spectro¬ 
scopic observations of IRAS 19135+3937, an SRd variable 
with an IR excess. It is often assumed that variability of 
such stars is due to pulsations, but there is no theoreti¬ 
cal study yet that would consistently reproduce all types of 
semi-regular variables. We showed that IRAS 19135+3937 
is a binary system with a circumbinary dusty disc. These two 
properties, taken together, can explain the observed light 
variations with a period of 127 days by an obscuration of 
the primary star by the inner edge of the disc during infe¬ 
rior conjunctions. Given the very smooth Kepler lightcurve, 
the photosphere of IRAS 19135+3937 is stable and does not 
show any pulsations. 

IRAS 19135+3937 presents a typical case of a class of 
pAGB binaries where circumstellar matter is confined to a 
disc: it is a low-gravity object, a single-line spectroscopic bi¬ 
nary, shows a depletion pattern of the refractory elements 
(albeit on a weaker side of the range), and has a jet-like 
outflow from the companion, based on the specific Ha pro¬ 
files. The commonality of the latter p henomenon was rea l- 
ized thanks to the HERMES survey ( Gorlova et al . 21113 !. 
Furthermore, we detected static emission lines of Fel and 
TiO that are most pronounced in the spectra taken near 
minimum light. We successfully fit the double-peaked Fe 
lines with a circumbinary Keplerian disc, that may be the 


first evidence of the star-disc interaction leading to the de¬ 
pletion pattern in pAGB atmospheres. The nature of the 
invisible companion remains unknown. We detected a sec¬ 
ondary set of spectral lines, but prove that it must be a 
reflected spectrum of the primary on the disc grains, rather 
than a spectrum of a physical companion. 

The single case of IRAS 19135+3937 presented here 
can not rule out pulsations for all SRd variables. What our 
study illustrates is that the presence of a circumbinary dusty 
disc may have profound effects on the light and RV curves, 
and therefore needs to be considered on a par with the more 
traditional factors, such as stellar eclipses and pulsations. 
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APPENDIX A: NEW MEASUREMENTS OF IRAS 19135+3937 


Table Al: UBV photometry in 2012-2013 


JD 

V 

B 

U 

to 


2456093.447 

11.730 

12.231 

12.521 

0.501 

0.290 

2456094.482 

11.716 

12.236 

12.485 

0.520 

0.249 

2456095.375 

11.731 

12.199 

12.522 

0.468 

0.323 

2456101.456 

11.674 

12.158 

12.454 

0.484 

0.296 

2456119.425 

11.299 

11.854 

12.194 

0.555 

0.340 

2456121.473 

11.204 

11.781 

12.147 

0.577 

0.366 

2456122.437 

11.180 

11.752 

12.086 

0.572 

0.334 

2456128.378 

10.955 

11.550 

11.940 

0.595 

0.390 

2456131.353 

10.859 

11.454 

11.874 

0.595 

0.420 

2456133.355 

10.781 

11.402 

11.854 

0.621 

0.452 

2456136.312 

10.719 

11.344 

11.773 

0.625 

0.429 

2456137.447 

10.702 

11.312 

11.758 

0.610 

0.446 

2456146.374 

10.615 

11.263 

11.732 

0.648 

0.469 

2456147.382 

10.628 

11.286 

11.743 

0.658 

0.457 

2456151.353 

10.644 

11.295 

11.765 

0.651 

0.470 

2456159.404 

10.793 

11.448 

11.895 

0.655 

0.447 

2456161.360 

10.844 

11.488 

11.951 

0.644 

0.463 

2456166.472 

10.996 

11.600 

12.012 

0.604 

0.412 

2456176.375 

11.318 

11.881 

12.234 

0.563 

0.353 

2456177.302 

11.372 

11.913 

12.269 

0.541 

0.356 

2456186.302 

11.597 

12.103 

12.365 

0.506 

0.262 

2456189.333 

11.598 

12.094 

12.382 

0.496 

0.288 

2456200.264 

11.734 

12.215 

12.504 

0.481 

0.289 

2456202.319 

11.740 

12.233 

12.545 

0.493 

0.312 

2456208.313 

11.722 

12.210 

12.487 

0.488 

0.277 

2456216.232 

11.687 

12.180 

12.504 

0.493 

0.324 

2456224.236 

11.695 

12.163 

12.481 

0.468 

0.318 

2456241.222 

11.352 

11.878 

12.192 

0.526 

0.314 

2456249.208 

11.017 

11.591 

11.964 

0.574 

0.373 

2456270.188 

10.613 

11.249 

11.660 

0.636 

0.411 

2456276.184 

10.674 

11.318 

11.760 

0.644 

0.442 

2456405.565 

10.669 

11.349 

11.831 

0.680 

0.482 

2456406.512 

10.682 

11.352 

11.852 

0.670 

0.500 

2456420.503 

10.908 

11.577 

12.045 

0.669 

0.468 

2456422.506 

10.945 

11.627 

12.059 

0.682 

0.432 

2456431.517 

11.229 

11.865 

12.216 

0.636 

0.351 

2456434.490 

11.315 

11.927 

12.278 

0.612 

0.351 

2456445.479 

11.641 

12.184 

12.455 

0.543 

0.271 

2456454.514 

11.743 

12.239 

12.577 

0.496 

0.338 

2456463.510 

11.760 

12.265 

12.546 

0.505 

0.281 

2456472.379 

11.748 

12.231 

12.528 

0.483 

0.297 

2456473.392 

11.725 

12.240 

12.547 

0.515 

0.307 

2456479.367 

11.676 

12.171 

12.510 

0.495 

0.339 

2456482.428 

11.609 

12.113 

12.434 

0.504 

0.321 

2456484.401 

11.568 

12.102 

12.409 

0.534 

0.307 

2456485.356 

11.558 

12.072 

12.402 

0.514 

0.330 

2456487.401 

11.504 

12.042 

12.324 

0.538 

0.282 

2456489.405 

11.437 

11.988 

12.332 

0.551 

0.344 

2456492.399 

11.343 

11.914 

12.309 

0.571 

0.395 

2456504.440 

10.860 

11.506 

11.933 

0.646 

0.427 

2456506.431 

10.801 

11.437 

11.891 

0.636 

0.454 

2456510.426 

10.679 

11.347 

11.801 

0.668 

0.454 

2456513.336 

10.649 

11.306 

11.773 

0.657 

0.467 

2456513.413 

10.637 

11.307 

11.737 

0.670 

0.430 

2456514.339 

10.627 

11.284 

11.787 

0.657 

0.503 


Continued on next page 
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Table A1 - Continued from previous page 


JD 

V 

B 

U 

to 

1 

U-B 

2456515.397 

10.606 

11.279 

11.744 

0.673 

0.465 

2456517.412 

10.594 

11.244 

11.751 

0.650 

0.507 

2456518.358 

10.592 

11.258 

11.755 

0.666 

0.497 

2456519.424 

10.572 

11.248 

11.729 

0.676 

0.481 

2456520.379 

10.575 

11.249 

11.697 

0.674 

0.448 

2456531.380 

10.635 

11.304 

11.776 

0.669 

0.472 

2456533.410 

10.670 

11.336 

11.819 

0.666 

0.483 

2456545.372 

10.965 

11.621 

12.043 

0.656 

0.422 

2456563.302 

11.615 

12.127 

12.454 

0.512 

0.327 

2456573.270 

11.749 

12.248 

12.519 

0.499 

0.271 

2456574.326 

11.738 

12.240 

12.508 

0.502 

0.268 

2456577.271 

11.781 

12.280 

12.602 

0.499 

0.322 

2456586.236 

11.807 

12.324 

12.649 

0.517 

0.325 

2456591.264 

11.790 

12.304 

12.603 

0.514 

0.299 

2456597.284 

11.731 

12.222 

12.575 

0.491 

0.353 

2456602.292 

11.702 

12.219 

12.505 

0.517 

0.286 

2456606.257 

11.687 

12.183 

12.477 

0.496 

0.294 

2456607.236 

11.695 

12.203 

12.510 

0.508 

0.307 

2456612.229 

11.688 

12.176 

12.463 

0.488 

0.287 


This paper has been typeset from a TJrjK/ file prepared by the author. 
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Table A2. Radial velocities in 2009-2013 


JD 

(d) 

Remain 
(km S -1 ) 

RVsec 
(km s _1 ) 

2454993.45906 

13.90 

-20.90 

2455003.50920 

3.29 

-22.24 

2455003.52830 

3.53 

-22.56 

2455003.54740 

3.62 

-19.50 

2455012.40789 

-6.29 

37.08 

2455012.42723 

-6.05 

33.09 

2455012.44660 

-6.23 

34.26 

2455025.56787 

-16.78 

25.22 

2455034.55947 

-18.74 

22.30 

2455061.44461 

-1.69 

N/A 

2455084.42500 

17.32 

-7.19 

2455423.61170 

-16.83 

14.74 

2455430.45750 

-12.33 

12.36 

2455501.33556 

13.52 

-21.61 

2455507.32911 

7.97 

-18.11 

2455650.71184 

-7.63 

33.15 

2455707.66632 

5.21 

-10.86 

2455714.50939 

6.58 

-13.75 

2455745.56824 

19.68 

-22.63 

2455761.54727 

7.06 

-24.93 

2455763.56835 

4.86 

-24.68 

2455773.61864 

-4.93 

30.83 

2455778.54334 

-8.67 

29.72 

2455784.49019 

-12.61 

29.06 

2455807.49337 

-17.13 

10.34 

2455828.42879 

0.22 

N/A 

2455837.44068 

-1.57 

N/A 

2455843.35922 

2.57 

-16.91 

2455873.34316 

20.68 

-17.48 

2455993.73419 

19.49 

-20.90 

2456010.70972 

13.30 

-20.14 

2456015.64342 

7.94 

-24.72 

2456033.67080 

-9.28 

27.91 

2456056.54496 

-16.41 

14.72 

2456082.59648 

0.01 

N/A 

2456087.52359 

-1.19 

N/A 

2456088.53941 

-1.90 

N/A 

2456103.40544 

11.66 

-14.23 

2456107.43456 

11.77 

-18.93 

2456117.48552 

21.29 

-18.95 

2456127.54549 

20.26 

-19.66 

2456136.41771 

14.06 

-22.11 

2456140.50515 

9.14 

-24.19 

2456144.61982 

4.28 

N/A 

2456146.48406 

3.85 

N/A 

2456152.45858 

-2.36 

30.74 

2456157.49963 

-7.44 

29.25 

2456165.48663 

-13.20 

27.11 

2456176.53346 

-14.96 

22.14 

2456178.58460 

-13.87 

21.02 

2456188.51469 

-9.01 

16.36 

2456194.44581 

-8.25 

N/A 

2456392.74729 

12.65 

-22.95 

2456441.69257 

-17.15 

11.11 

2456457.45702 

1.96 

N/A 

2456468.48879 

2.90 

-14.97 

2456477.43936 

2.73 

-15.64 

2456488.56607 

17.36 

-15.94 

2456512.62342 

18.19 

-22.41 

2456544.46648 

-11.58 

26.57 

2456566.49203 

-11.40 

16.88 


HERMES radial velocities for the main and secondary components in the cross-correlation function of IRAS 19135+3937, as described 
in Sec. 1 5.1 1 “N/A” stands for ’not available’ and designates phases where it was not possible to disentangle the secondary component. 





